Entanglements in Systems with Multiple Degrees of Freedom 
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In this letter we present the entanglement properties of the spin-orbital coupling systems with mul- 
tiple degrees of freedom. After constructing the maximally entangled spin-orbital basis of bipartite, 
we find that the quantum entanglement length in the noninteracting itinerant Fermion system with 
spin and orbit is considerably larger than that in the system with only spin. In the SU(2)<g>SU(2) 
spin-orbital interacting system, the entanglement, expressed in terms of the spin-orbital correlation 
functions, clearly manifests the close relationship with the quantum phases in strongly correlated 
systems; and the entanglement phase diagram of the finite-size systems is in agreement with the 
magnetic and orbital phase diagram of the infinite systems. The application of the present theory 
on nucleon systems is suggested. 
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The total spin wavcfunction of two entangled particles 
with spin-1/2 can not be expressed as the direct product 
of the wavefunctions of the two individual particles, and 
the measurement to the spin of one individual particle in- 
evitably brings out the information of the spin states of 
the other particle. This property provides wide potential 
application in quantum computer and quantum telepota- 
tion. Entanglement properties in quantum many-particle 
systems have received great interests in recent years since 
it is realized that, on the one hand, the entanglement may 
characterize the quantum correlations between particles 
0>0> or even the quantum phase transitions in some sim- 
ple models (3, 0, H El ; on the other hand, the realization 
of quantum computation in solid devices also naturally 
raises the question, e.g. how these electrons entangle 
with each other in the presence of strong correlations. 
Up to date, most of studies are concentrated on the spin 
entanglement of many-particle systems. While in quan- 
tum many-particle systems, besides the spin freedom de- 
gree, the particle may possess other degrees of freedom 
(DOF), such as the orbit, the position and the momen- 
tum, etc. Especially in the strongly correlated electronic 
systems, the orbital DOF is important for the relatively 
localized 3d or 4f electrons jj. Little was known about 
the entanglement properties of the particles with mul- 
tiple freedom degrees. In this letter we present the en- 
tangled properties of strongly correlated particles with 
multiple DOF. After constructing the maximally entan- 
gled states, the Bell basis, of two particles with spin and 
orbit, we first demonstrate that the quantum entangle- 
ment length in itinerant electron system with two DOF 
becomes significantly larger than that in the system with 
only spin. In the strongly interacting spin-orbital systems 
with STJ(2)(g)SU(2) symmetry, the entanglement evolves 
with the quantum phases, and appears abrupt changes 
near the critical points of the quantum phase transitions. 
The entanglement phase diagram highly coincides with 
the magnetic and orbital phase diagram. Finally we dis- 



cuss the possible application of the present theory on the 
nucleon systems with multiple DOF. 

Considering the quarter-filling (v—l/A) quantum 
many-particle systems with spin S=l/2 and twofold- 
degenerate orbital DOF, we introduce the pseu- 
dospin operator f to describe the two orbital states: 
T=^2 ab C\(j a bCb, here a is the Pauli matrix, and the or- 
bital indices a, b run over 1 and 2. As well known, an en- 
tangled state is unfactorizable, and a bipartite with max- 
imally entangled state is completely indivisibility. For a 
bipartite with two electrons at site A and site B with spin 
DOF, the Bell basis consists of four states \ip{ 2 ) ab — 

(| U) ± I IT)) /y/2, and \^ a )ab = (| TT) ± I U» fy/2. 
These Bell states with maximized entanglement are the 
foundation of quantum information theory and quantum 
computation [2j . In the presence of both spin and orbital 
DOF, we find that the possible Bell basis is separated 
into two independent groups. 

(i) The first group is composed of the direct product of 
the entangled spin and the entangled orbital parts, and 
each part is a Bell sub-basis, i.e., maximally entangled. 
The 16 basis wavefunctions in the present situation can 
be readily constructed: \i/j st )ab = \4> s ) BeU <g> |^ T ) Be ", 
where \ip s } Bel1 and \ip T ) Bel1 represent the spin and orbital 
Bell states, respectively, e.g. 

\r T h,2 = (|TT>±IU»®(|n>±|22» (i) 

etc. Similar to the spin case, the 16 states in this 
group are the common eigenstates of the operator set 

^A^B T A T B^ VaOb t A t Bi a A°B T A r B^ °A°B T A r B)- 0nC 

can use these operations to generate the 16 states in this 
group without difficulty. 

(ii) The second group of the Bell basis is consisted of 
16 spin-orbital states which arc the linear combination of 
independent spin-orbital parts, e.g. 

|<ni,2 = (|H) T <8|^)±|22) T (»|^))/V2 (2) 
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etc., here is the i-th spin Bell subbasis. We no- 
tice that as the candidates for the maximally entangled 
basis, each group of these basis is orthonormalized and 
complete set. 

In order to justify if these are the maximally entangled 
states, the degree of entanglement of these states is 
quantified. For a many-state bipartite, the negativity is a 
good operational measurement to quantify the entangle- 
ment since it is monotonic under the local operation and 
the classical communication, and it vanishes if two sub- 
systems are not entangled ^a,b — (\\p Ta \ \ — l) /2, 
where ||/0 Ta || is the trace norm of the partial transpose 
of the density matrix of subsystem A versus to the B 
and satisfies: ||/3 Ta || = Try/ p TA (p TA )+ . Interestingly, 
since spin and orbit are two independent DOF, we could 
introduce the spin and the orbital sub-entanglements 
to quantify the entanglement degree of each DOF after 
trace out the other DOF: K S M = - l)/2, where 

PsM = Tr T { s \(p). We find both of the two groups of the 
basis have the maximal entanglement with the von Neu- 
mann entropy Eab — 1 an d the negativity ^ab — 1-5. 
However the spin and the orbital sub-entanglements of 
these two subgroups are different: in the first group, 
the spin or orbital maximal sub-entanglement of each 



state is N 



AB — 



0.5; whileas in the second group, the 
spin or the orbital sub-entanglement of each state 
completely vanishes. Such significant difference in 
the sub- entanglement reflects the distinct character of 
these two group basis: the former is composed of the 
independent spin or orbital Bell states, while the latter 
of spin-orbital indivisible bipartite. 

A. Entanglement in Itinerant Fermion Systems 
Before investigating the entanglement in the spin-orbital 
interacting system, we first explore the bipartite entan- 
glement in non-interacting itinerant electron system with 
spin and orbital DOF. The spin entanglement of two elec- 
trons in noninteracting Fermi gas with single spin DOF 
was known |lfl 111! Il2l | . In the system with only spin, 
the entanglement vanishes when the spatial separation 
of two electrons is larger than a characteristic length r°, 
r~ « \.&/k F , where k F is the zero-temperature Fermi 
momentum. In the spin-orbital itinerant electron sys- 
tem, we find that the spin-orbital entanglement length of 
two electrons becomes considerably large, and the sub- 
entanglement of the bipartite for individual DOF is zero 
even in the most entangled spin-orbital states. 

At zero temperature in the ground state, the elec- 
trons occupy all the levels below the Fermi surface 
IV'o) = rL,r,|fc|<fc F c k.s.rl0)> where the Fermi wavevec- 
tor k F satisfies k F = (3n 2 N/2V) 1/3 . From the den- 
sity matrix p~N — |V , o)(V'o|j the two-particle density 
matrix p\i is obtained through tracing over other par- 
ticles pi 2 = Tr{ PN (f>+(2')cj)+(V)(j)(l)(j){2)} [H, where 
1 = (ri,<7i,Ti), ri is the position vector, and the oper- 



ator (f)(1) = l/VV"^2k e ~* k ric k,s,T- Taking into account 
the spatial diagonal element, we obtain the reduced two- 
particle spin-orbital density matrix, 

/[16-4/(r) 2 ] (3) 

where f(r) = 3(sinx — xcosa;)/x 3 with x = k F r. For 
this mixed-state density matrix, the negativity is larger 
than zero as f(r) > 1/2, 



/(rf 



f{rf 



(4) 



The negativity of two electrons as the function of the 
distance of k F r at T=0 K for the spin-orbital and spin 
systems is shown in Fig. 1. We find that in the present 
spin-orbital system the quantum entanglement length r e 
is about 2.4/fcp, considerably larger than that in the spin 
system. Considering the shrinking of the Fermi energy in 
the quarter-filling spin-orbital system, the entanglement 
length is about 1.7 times larger than that in the spin sys- 
tem. In the present itinerant electron system, the two- 
particle spin-orbital density matrix (3) can be expressed 
as: p = [4(1 - /(r) 2 )//16 + 3/(r) 2 p']/(4 - /(r) 2 ), where 
p' is the density matrix composing of all spin-orbital two- 
particle antisymmetric states and I is the 16-order unit 
matrix. As a contrast, the 4-order density matrix in the 
spin system is Po = [(2(1 - /(r) 2 )/ /4 + f(r) 2 p' ]/(2 - 
/(r) 2 ), where p' is the density matrix of spin singlet 
state and Iq is the 4-component unit matrix. Evidently, 
the unit matrix represents completely mixed disentan- 
gled state, so the more the weight of the unit matrix in p 
is, the more disentangled the state is. At sufficient large 
distance k F r, the weight of the unit matrix, or the dis- 
entangled component, in the spin-orbital density matrix 
p is much smaller than that in the spin density matrix 
po, leading to a significantly large quantum entanglement 
length in the spin-orbital system. 

For two particles located at the same position, their 
spatial wavefunctions are the same. Due to Pauli princi- 
ple, the spin and orbital wavefunctions must be antisym- 
metric, and it does not contributes to the unit matrix. 
Thus the entanglement of the two electrons at is maxi- 
mal at r = 0, as shown in Fig.l. Meanwhile, the spin- 
orbital wavefunction of the two electrons is composed of 
six antisymmetric spin-orbital substates, e.g. | ^1}®\ip~), 
etc., here |Vv ) i s the orbital singlet state. These compo- 
nents are not the Bell basis wavefunctions, so p, though 
is equal to p' , is still a mixed-state density matrix, and 
the negativity is much smaller than the maximal value 
1.5. Comparing to the spin system, the entanglement 
length and the Fermi velocity in the present system sat- 
isfy r e rs 1.7?"g and v F — v F /2s , respectively. Therefore, 
the coherent time of the entangled electrons t c ~ r e /v F 
is about 2.1 times larger than that of the spin system. 
This property is also expected to valid for the localized 
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electron systems. Consequently, the present system with 
spin and orbit DOF is more favorable than the system 
with only spin in the realization of quantum computation 
and quantum telepotation . 
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FIG. 1: Negativity as the function of spatial distance Uft 
of two electrons. N and JV° represent the negativity of two 
electrons in the spin-orbital system and in the spin system, 
respectively. 



It is interesting to explore the variation of the spin sub- 
entanglement of two electrons with the spatial distance. 
The reduced spin density matrix becomes, 



1 



Pl2 



f(r) 



al (4-2/(r)')f 



/(rfU-XVa (5) 



after tracing over the orbital DOF, where \4>j) is the spin 
singlet state. We find that the spin-orbit density matrix 
of the two electrons at the same position, /of 2 ( r = 0), 
is equal to that in the spin system at the entangled 
distance, p°(r e ), indicating that even at the maximally 
entangled spin-orbital states, the spin of two electrons in 
itinerant electron system is disentangled. The vanishing 
spin sub-entanglement in itinerant electron systems 
arises from the enlargement of Hilbert state space, which 
greatly reduces probability of forming spin singlet state. 
The spin correlation functions, (SiSj), also manifest the 
disentanglement of the spin DOF. We find at r = 0, 
(SiSj) = —1/4, showing that the two electrons are 
classical AFM correlated. With the increase of the 
separation between two electrons, the classical AFM 
correlation becomes smaller and smaller. Therefore the 
spin sub-entanglement in any spin-orbital state vanishes 
at any distance in itinerant electron systems. 

B. Entanglement and Phase Diagram in Localized Elec- 
tron Systems 

To explore the relationship between the variation of the 
entanglement and quantum phase transitions in many- 
electron systems, we study the entanglement in the lo- 
calized spin-orbital interacting systems in what follows. 
Consider a one-dimensional system with strong spin- 



orbital correlation 0,0]. Its Hamiltonian reads: 

H = /][Sj • Sj + lTj • T l+ l + JgSj • Sj + l + J T Ti ■ T l+ l 

i 

+B s Sf + B T T*] (6) 

where J s and J T are the spin and the orbital exchange 
constants, respectively, and B s (B T ) is a small mag- 
netic (orbital) field. In the absence of the small exter- 
nal fields, both the spin and the orbital parts in Eq.(6) 
satisfy the SU(2) symmetry, the so-called SU(2)<g>SU(2) 
model. The periodic condition on Eq.(6) implies that 
the entanglement between arbitrary two sites is an uni- 
form function of distance. The Hamiltonian commutes 
with the z-components of the total spin and the total 
orbital operators, [H, S z ] = and [H, t z ] = 0, where 
S z = J2i &i an d tZ — J2i T i- S° each eigenstate of the 
Hamiltonian is also the eigenstate of S z and r 2 . One 
can use the index (S z ,t z ) to characterize the main four 
phases of this model in different regions of the exchange 
constants J s and J T : Phase I with (S z =N/2, r z =N/2), 
Phase II with (N/2,0), Phase III with (0,N/2) and Phase 
IV with (0,0), as shown in Fig. 2. In the present situation 
the spin-orbital entanglement can be explicitly expressed 
in terms of the spin-orbital correlation functions. Using 
the reduced density matrix of two nearest-neighbor parti- 
cles, one can analytically express the negativity H[j as 
the nearest-neighbor spin-spin, orbital-orbital and spin- 
orbital correlation functions in each phase. For instance 
in the phase IV, the negativity reads: 



[*,»+!] 



2\d\ 
+9- 



|c-ft| + |c + /i| + |6-/| + |6 + /| 
— g — e + k\ + \a — g + e — fc| + \a 
i-k\ + \a + g + e + k\-l}. (7) 



with 



a/b = 


((\ + s:s: +1 )(±±t z t z +1 )) 




c/d = 


((J-S?Sf +1 )(i±7?7? +1 )> 




e/f = 


((±±T*T* +1 )(S*Sf +1 +S?SV_ 




9/h = 


(( 1 -±s:s: +1 MT: +1 +TyT? + 


l)> 


k = 


It QX qX j_ qV qV \( X X , 


V 
T-T 

i i 



i+ i)>. (8) 

Obviously, in the fully polarized orbital phase III, the 
negativity has a simple form: = —1/4— (S;Si +1 ). 

Similarly we obtain the negativity in the completely po- 
larized spin phase II, = —1/4 — (riri+i). The 
negativity vanishes in the fully polarized spin and orbital 
phase I. Therefore, the direct relation between the quan- 
tum entanglement and the spin and orbital correlations 
in strongly correlated electronic systems is thus estab- 
lished, and it may provide many interesting information 
hidden in strongly correlated systems. 



4 



Thus one could clearly find that the variation of the 
entanglement closely relates to the transition of quan- 
tum phases in the present system: at the critical point 
of the quantum phase transition, the entanglement ex- 
hibits a discontinuous change. The entanglement phase 
diagram of a 4-site finite system is shown in Fig. 2. In 
the phase diagram Fig. 2a, we adopt and the sub- 

s(t) 

negativity \^ to characterize these different quan- 
tum phases. Very interesting, we find that the finitc- 




-1 -0.5 0.5 1 -0.1 0.2 0.5 0.8 

FIG. 2: Entanglement phase diagram of 4 particles in (a) and 
dependence of negativity on magnetic exchange constant J s 
at J T = 0.25 and J T = 1.0 in (b). Phases I, II, III and IV 
represent the FM and FO state with (S z =2, r 2 =2), the FM 
and AFO state with (2,0), the AFM and FO state with (0, 2), 
and the AFM and AFO state with (0, 0), respectively. Phase 
IV differs from IV' in the spin-orbital correlation function, see 
the address in the text. 



changes of entanglement of two electrons near the phase 
boundary are in accordance with the critical points of 
the quantum phase transitions, as clearly shown in Fig 
2(b). At 3 S —3 T =1/A and B s =B T =0, the system exhibits 
SU(4) symmetry |17|. then the ground state is highly 
degenerate and the negativity comes to the maximum. 
These results clearly demonstrate that the entanglement 
is closely related with the quantum phase transitions in 
strongly correlated systems. It is worthy of pointing out 
that we introduce the small magnetic and orbital fields 
to lift the degeneracy of Phases II and III. In the ab- 
sence of these small fields the entanglement is uncertain, 
as discussed by Qian, et al (lif. 

It is interesting that the present theory is also appli- 
cable for the nuclear systems, in which the nucleons usu- 
ally possess more than one DOF, such as spin, isospin, 
etc. According to our results for the itinerant Fermion 
systems, the entanglement length of the nucleons is sig- 
nificantly larger than one expects. We anticipate more 
interesting entanglement properties in the systems with 
multiple degrees of freedom will be uncovered in the fur- 
ther studies. 

Authors thanks for the useful discussion with Z. S. Ma. 
Supports from the NSF of China and the BaiRen Project 
from the Chinese Academy of Sciences (CAS) are appre- 
ciated. Numerical work was performed in CCS, HFCAS. 



site entanglement phase diagram in Fig. 2(a) is almost 
identical to the magnetic and orbital phase diagram in 
the one-dimensional infinite lattice [14L Il5| . indicating 
that the quantum information of infinite systems mani- 
fests through the entanglement of proper finite clusters. 
The phase diagram in Fig. 2a can also be characterized 
by the spin-orbital index (S z , t z ). Among these phases 
in Fig. 2a, the states in Phase I are disentangled, corre- 
sponding to fully polarized spin and orbital states with 
(S z —2, t z —2), i.e. the ferromagnetic (FM) and ferro- 
orbital (FO) state. In Phase II, the negativity satisfies 
N[i,i+i]=N[i and H^ i+1 j=0, indicating the states in 
Phase II are FM and antiferro-orbital (AFO) with (2,0). 
Similarly, Phase III corresponds to the AFM and FO 
states with (0,2). Both Phases IV and IV correspond 
to AFM and AFO with (0,0), while these two phases 
are distinguished by the spin-orbit correlation functions. 
In Phase IV, (Si • S^+ir^ • Ti + i)<0, and in Phase IV, 
(Si • Si + iTi ■ -r; + i)>0. This difference arises from the fact 
that when either J s or J T becomes large while the other 
one is small, (Si ■ Si+i) < and (t^ ■ r^+i) < 0. Due 
to the quantum fluctuation, (Si ■ Si + iTi • Ti + i) < 0, in 
Phase IV. When both J s and J T become large enough, 
the strong AFM and AFO correlations polarize the spins 
and orbits, and lead the spin-orbital correlation function 
(S,; • Si + iTi • Ti + i) to transition from negative to positive, 
as we see the phase IV in Fig. 2a. The discontinuous 
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